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The transmission poles of N number of identical Dirac delta potentials placed periodically in 
one-dimension are examined in the complex-energy plane. The numerical results show that the 
imaginary part of the poles scales with l/N. An approximate form of the poles is derived which 
supports the scaling behaviour of the poles found numerically. It is shown that the imaginary part 
of the poles are proportional to their real part for the poles close to the ends of the bands. 



I. INTRODUCTION 

In recent years, quantum transport in mesoscopic systems has become a rapidly developing research are In the 
pioneering papers by Landauer and Biittiker it was shown that the conductance of such devices can be related to the 
transmission coefficients depending on the scattering properties of the systemEm. In electron waveguide structures the 
resonant behaviour and the transmission phenomena have become an extensively studied problemD'a. The transmission 
resonance caused by the quasi-bound states are in close connection with the poles of the transmission amplitude in 
the complex-energy plane. The quasi-bound state at energy Eq and with decay time r = fi/T corresponds to a simple 
pole in the transmission amplitude t(z) at the complex-energy z = Eq — iT/2 and the energy dependence of the 
transmission probability T(E) — |i(-E)| 2 for r <C Eq is given by the well-known Breit-Wigner formulatJo 

r 2 /4 

T(E) = 4 . (1) 

(£-£ ) +r 2 /4 

The pole structure of the transmission of a one-dimensional periodic array of Dirac delta potential (Dirac comb) is 
investigated in this paper. The transmission amplitude of N unit cells is expressed by the transmission of one unit 
cell. The poles of the transmission amplitude are calculated numerically and a definite structure of the poles in the 
complex energy plane shows up. An analytical expression is derived for the imaginary part of the poles and it is 
shown that it scales with A -1 , where N is the number of unit cells in the system. The agreement with the numerical 
results are excellent. 

The rest of the paper is organised as follows. In section [n| using the transfer matrix method the transmission 



coefficient of N number of unit cells is derived for the case of Dirac comb. In section III the numerical calculation 
of the poles are presented. In sectio n [IV| an approximate expression of the poles are derived and compared with the 
numerical results. Finally, in section M our conclusions are given. 



II. THE TRANSFER MATRIX APPROACH 



In this section the transmission of an one-dimensional periodic array of ./V Dirac delta potentials is considered. In 
one-dimension the single-channel description of the scattering process is valid. In this case Sprung et. al. derived 
the following compact expression for the transfer matrix of N number of periodically placed identical cells in terms 
of the 2x2 transfer matrix M of a single-cell: 

M N = — - — [Msin7V<I> — Ism (TV — 1)$], (2) 
sin $ 

Similarly, the following simple expression holds for0 the transmission probability Tjy = |tjv| 

1 _ sin 2 iV$ / 1 
]t^^ + sin 2 $ \W 

where t\ is the transmission amplitude of a single-cell and the Bloch phase $ associated with the infinite periodic 
potential is given by 



where cos $ = 1 /2 Tr M. 
of the A-cell array: 
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cos$ = Re (1/ti), 



(4) 



Note that the transmission amplitude depends on the energy of the electron involved in the scattering process. 

These results can be applied for a one-dimensional periodic array of N Dirac-delta potentials. The transmission 
amplitude t\ for the case of unit cell containing a Dirac delta potential of strength A and located in the centre of the 
unit cell is given by 

* = rrW" < 5 » 

where k = y/E, the length of the unit cell is denoted by d and E is the energy of the electron. Hereafter we use the 
h = 2m = 1 units. Using Eq. (^) we find 

cos$(E) = cos fed + sin fed. (6) 
2k 

Tha quasi-bound states of the system are given by the poles of the transmission amplitude in the complex energy 
planet]. Hence, for the TV-cell system l/T N = which, using Eq. (ft), can be written as 



sin 2 iV$ / 1 



sin 2 * VlTi 



0. (7) 



Thus, substituting Eq. (|^) into the above equation one finds that for Dirac comb the complex energy E of the 
quasi-bound states is the solution of the following equation: 

AE _ sin 2 N$(E) 
~ A^" ~~ sin 2 ( } 

where <fr(E) is given by Eq. (ph. 



III. NUMERICAL RESULTS 



We now present the numerical solution of Eq. (^) . In Fig. |l| the poles in the complex energy plane are shown for 
different number of Dirac delta potentials. 
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FIG. 1. The poles of the transmission amplitude of the Dirac comb for N = 8, 16, 32 number of Dirac deltas (the poles are 
denoted by plus signs, crosses and stars, respectively). The strength of the Dirac delta potentials, A = 10.0, and the distance 
between them, d — 4. 
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It can be seen that increasing the number of Dirac delta potentials, A in the array the poles move towards the real 
energy axis. It is also clear from the figure that the real part of the poles are separated into allowed and forbidden 
bands. The formation of the band structure of the infinite periodic array (Kronig-Penney modeH) can be seen even 
for A = 8. In each allowed band there are A — 1 number of poles. The width of the forbidden bands decreases with 
increasing the real part of the complex energy. The real part of the poles, E™ can be given by 

$(£f ) = lir/N, (9) 

where I = 1, ... A — 1 and $(E) is given by Eq. (^|). Hence, the real part of the Ith poles in the mth band, E" 1 is 
determined by 

cos ki m d + ^ sinki m d — cos(lir/N), (10) 

where fe; m = JEV 1 . 



IV. SCALING OF THE POLES 



In order to see the scaling behaviour of the poles an approximate solution of Eq. (g) for the quasi-bound states is 
derived. 

The solution of Eq. (||) is a complex energy in the form of E = E™ — iY7 l /2, where E™ is the A — 1 number of 
distinct solutions of Eq. ([Io|) indexed by I in the mth allowed band and the decay time of the quasi-bound state is 
h/Tf 1 . Assuming TV l /2 <C E™, the Bloch phase $(£') defined by Eq. (||), can be expanded in first order in T™ around 
EV l : 



where 



Finally, Eq. (0) becomes 



$(£)«$(£ )+<5$(i?r\rn, 

. d$(E) 



S^iE™^™) 



dE 



IT 



4(E; m + «T|"/2) _ sin 2 N(^>(E'[ n ) + id^Ef 1 , T™)) 



A 2 



sin 2 ($(^ m ) + i5<I>(E™, TV 1 )) 



(11) 



(12) 



(13) 



Substituting Eq. (Q) into Eq. (g_3|) and neglecting 5$ in the denominator of the right hand side of Eq. ( |l3|) results in 
an implicit equation for T™ of the quasi-bound state: 



pm 1 

= — arsmh 

2 N 



2y/E™ + i-L-sia(!K/N) \ ( d $(E) 



dE 



(14) 



where the last factor can be determined from Eq. 

d$(E) d 



dE 



4k sin $ 



2 H ] sin kd - 

Ed , 



A cos kd 



and k = yE. For TV 1 < Ef 1 the implicit Eq. (|14j) can be solved iteratively for T\ n . Substituting T[' 
right hand side of Eq. (n~3) yields 



rr 



l 



A (A/fc 2 



4fc im sin (/tt/A) 



+ 2c?) sinfc; m c? - {Xd/kim) coski m d 



irsinh 



2k lm sin (ln/N) 
A 



(15) 
into the 

(16) 



The approximate solution of Eq. (|J) as the complex energy, E™ — iT™ /2 of the Ith. quasi-bound state in the mth band 
is given by Eqs. @ and (|l|). 
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In Fig. g the numerically determined poles are shown for TV = 8, 16, 32 number of Dirac delta potentials (this is 
an enlarged portion of Fig , [j]) in the range of 66 < Real{E} < 75 together with the approximate analytical solution 
given by Eqs. ( |l0| ) and (|l^) when Ef 1 is taken as a continuous variable. 




FIG. 2. One band of the poles (an enlarged portion of Fig. [I] for 66 < Real{E} < 75) and the plot of Eq. ([TJ|) are shown. 
The number of Dirac deltas are N — 8, 16, 32 (the poles are denoted by plus signs, crosses and stars, respectively). 

One can see that the envelope of the poles is well approximated by Eq. (|l6|). The agreement between the analytic 
and the numerical solutions is excellent. Similar agreements are found in the other bands. The approximate solution 
becomes even better for increasing N. A small deviation appears only in the middle of the band but increasing N it 
is decreasing. From Eq. ( |l6| ) one can see that the imaginary part of the complex energy of the poles are proportional 
to 1/N. It is easy to show that the factor multiplied by 1/N in Eq. (|l6| ) does not depend on N in the leading order 
of 1/N. Therefore, for large N the imaginary part of the poles scales as 



1/N 



(17) 



for all the poles. 

Using the numerical solutions of Eq. (||) in Fig. || AT" l /2 is plotted as a function of the real part of the poles for 
N = 8, 16,64. As it is seen from the figure the numerically found poles follows the scaling behaviour of the poles 
given by Eq. (|l7|). Again a very small deviation arises only at the middle of the bands for N < 16. 




FIG. 3. NT™ /2 as a function of Real{E} for the case of N = 8, 16, 32 number of Dirac deltas. The same parameters of the 
system and signs are used as in Fig. |l|. 

Using arsinh(x) k x for i < 1 (in this case sin $(_£"') <C 1) in Eq. (|l6|) one can show that 



rr 



e. 



8 



sin 2 (lir/N) 



1 XN {X/kf m + 2d) sin k lm d - (Xd/ k lm ) cos k lm d ' 



(18) 
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for fixed N and I <C N or TV — I <C iV (these poles are close to the end of the band m). Again it can be shown that 
the factor multiplied by E™ in Eq. ( |l8|) is weakly dependent on the band index m for those values of E™ when all 
the terms in the denominator of the right hand side of Eq. (|l^) is negligible except for the term 2dsmki m d. Thus, 
in the complex energy plane the poles at the end of each band move down from the real axis as the band index m 
increases. This tendency can also be seen clearly from Fig. ^. However, from Eq. ( |l8| ) a more rigorous statement can 
be made, namely at fixed N and poles with I close to the end of the band Tf 1 ~ EJ n . In Fig. || the imaginary part of 
the poles, T™ calculated numerically is plotted as a function of the real part of the poles, -E™ for I = 1 as well as the 
analytic form given by Eq. JI§|). 
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FIG. 4. r^/2 as a function of -E™ for N = 32 and I = 1, and the plot of Eq. (|l^). The same parameters of the system are 
used as in Fig. [j]. 

The numerically calculated poles are in good agreement with the analytic form given by Eq. (|l8| ) even for TV = 8. 
The poles move down from the real axis in a way that their imaginary part are proportional to their real part for 
those poles which are close to the ends of the bands. 

V. CONCLUSION 

We have studied transmission properties of the one-dimensional periodic array of Dirac delta potential. Using the 
transfer matrix method a closed form of the transmission amplitude is given for finite number of Dirac delta potentials 
(finite number of unit cells in the Dirac comb). The quasi-bound states, which affect on the transmission properties 
of the system, are related to the poles of the transmission amplitude. The pole structure of the system is investigated 
for different number of unit cells in the Dirac comb. 

The numerical solutions reveal a definite structure of the poles. The poles are well separated into allowed and 
forbidden bands much in the same way as in the Kronig-Penncy model but the 'pole-spectrum' is not continuous. 
Each band contains TV — 1 number of poles, where N is the number of unit cells in the system. To check our numerical 
results an analytical form for the imaginary part of the poles has been derived. Although this is an approximate 
expression the comparison with the numerical results shows an excellent agreement especially for large N. From our 
analytic solution it has also been proved that the imaginary part of the poles is proportional to X/N. Furthermore, 
it has been demonstrated that for the poles close to the ends of the bands the imaginary part of the poles are again 
proportional to the real part of the poles. 
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